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Abstract. We prove that the double layer potential operator and the gradient 
of the single layer potential operator are L2 bounded for general second order 
divergence form systems. As compared to earlier results, our proof is independent 
of well posedness for the Dirichlet problem. More precisely, we use functional 
calculus of differential operators with non-smooth coefficients to represent the 
layer potential operators as abstract bounded Hilbert space operators. In the 
presence of Moser local bounds, in particular for real scalar equations and systems 
that are small perturbations of real scalar equations, these abstract operators are 
shown to be the usual singular integrals. Our proof gives a new construction of 
fundamental solutions to divergence form systems, valid also in dimension 2. In 
general, the layer potential operators are shown to depend holomorphically on the 
coefficient matrix A € Loo, showing uniqueness of the extension of the operators 
beyond singular integrals. 



1. Introduction 

This paper concerns the classical boundary value problems for divergence form 
second order elliptic systems 

m. 

^ divA'^Vu^ = 0, i = 1, . . . , m, 

i=i 

for a vector valued functions u = (u^)^^-^ on the upper half space R^"^" := {(t,x) G 
R X R" ; t > 0}, n,m > 1, with boundary data in L2(R"). In general, we 
only assume that the coefficients A = (A*-')™-^]^ are uniformly bounded and accre- 
tive. Unless otherwise stated, we assume that A{t, x) = A{x) is independent of the 
transversal direction t. (Accretivity, or more precisely strict accretivity, is defined 
in (jl]) below.) However, we do not assume that A is real or symmetric. 

By scalar coefficients, or equation, we mean that A'^^ = for z 7^ j. For technical 
reasons we consider systems where the functions are complex-valued, and thus 
A'^^[t,x) G C{C^~^^). However, working at the level of systems of equations of arbi- 
trary size, complex coefficients are no more general than real coefficients. Indeed, 
using the relation C = R^ we see that any system of equations with complex coeffi- 
cients of size m can be viewed as a system of equations with real coefficients of size 
2m. 
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For an 1 + n-dimensional vector /, we let denote the normal/vertical part 
(identified with the corresponding scalar component), and write /y for the tan- 
gential/horizontal part. Similarly, we write Vy, divy and curly for the differential 
operators acting only in the tangential/horizontal variable x. To ease notation, we 
use the Einstein summation convention throughout this paper. Sometimes we shall 
even suppress indices 

A classical method for solving the Dirichlet problem is to solve the associated 
double layer potential equation on the boundary R"^. In our framework, the method 
is the following. Let T(^t,x) = {'^^('t x))Tj=i be the fundamental solution for div74*V 

in Ri+" with pole at (t,x), that is div(AJ*)*Vrff . = J ^ " ^' ^nd let 

^•^A*^^{tx) ■~ ii^''^)*^^'{t x)) ^ denote its (inward) conormal derivative. 
Given a function h : R" — )■ on the boundary, define the function 

Wix) := / {-d,^^T^^JO,y),h^iy))dy, (t,x) G R^+", 

where —di,^, is the outward conormal derivative. The function u{t,x) := Vth{x) 
then solves the equation divAVw = in R+'^"', and has boundary trace 

Vh^x) :=\\mJ^^{-d,^S';,^){0.y),h^{y))dy. 

Finding the solution u with Dirichlet data (p : R" — t- C"* on the boundary, then 
amounts to solving the double layer equation 

Vh = LP 

for h, which then gives the solution u{t,x) = Vth{x). In the case of smooth coef- 
ficients A, it is well known that the operator V is well defined and is |/ plus an 
integral operator. For general systems with non-smooth coefficients, as considered 
in this paper, the double layer potential operator V is beyond the scope of singular 
integral theory. 

Similarly, the single layer potential is used to solve the Neumann problem. See 
Section [71 In this introduction, we focus on the double layer potential and the 
Dirichlet problem. 

During the last years, new results on boundary value problems for more general 
non-smooth divergence form systems have been proved. In particular, there has been 
two seemingly different developments, one based on singular integrals (S) and one 
based on functional calculus (F). The purpose of this paper is to demonstrate that 
the singular integral operators used in (S) actually are special cases of the abstract 
operators used in (F). 

(S) In the paper by Alfonseca, Auscher, Axelsson, Hofmann and Kim, it was 
proved in particular that boundedness and invertibility of the layer potential 
operators for coefficients Aq implies boundedness and invertibility of the 
layer potential operators for coefficients A whenever \\A — Ao||oo is small, 
depending on Aq. Here Aq and A are assumed to be scalar and complex, 
and such that De Giorgi-Nash local Holder estimates hold for solutions to 
these equations. Boundedness here includes square function estimates. This 
boundedness and invertibility result was shown to hold for real symmetric 
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coefficients, and the result was also known for coefficients of block form and 
for constant coefficients. 

During the writing of this paper, Hofmann, Kenig, Mayboroda and Pipher [U] 
have proved Lp well posedness, for some p < oo depending on A, of the 
Dirichlet problem for general scalar equations with real and t-independent 
coefficients. From this they deduce, in [HI Cor. 1.25], boundedness in L2 (but 
not invertibility) of the layer potentials for general real scalar equations and 
small complex perturbations of such, by inspection of the proofs in [1]. 
(F) Auscher, Axelsson and Mcintosh [3] proved that the L2 Dirichlet (and Neu- 
mann) problem is well posed for systems with coefficients A which are small 
Loo perturbations of Hermitean, constant or block form coefficients. Instead 
of the double layer potential operator V above, this used an operator T) on 
L2(R") defined by functional calculus from an underlying differential oper- 
ator on R". More precisely, this used a self-adjoint ffist order differential 
operator D and a transformed multiplication operator B formed from the 
coefficients A, to construct a solution 

u\t,x) = Vth\x) ■= {btiBD)h)\{x), (t,x) G R^+", 

{e~*^, Rez > 0, 
' ' is applied to the operator BD. 

0, Re2; < 0, 

Here BD, and therefore bt{BD), acts on C™*-^"*""''-valued functions h on R". 

Both works [H [3] build on harmonic analysis developed for the solution of the Kato 
square root problem by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian 
However, the approach (F) is more general. On the one hand (S) uses De Giorgi- 
Nash local Holder estimates, which holds for real scalar equations, and small L^o 
perturbations of such, but not for general A. On the other hand, (F) proves that 
in fact is L2-bounded for any t-independent and uniformly bounded and accretive 
coefficients A; it is only invertibility oiT> := limi_^o+ which may fail. 

Unlike V, the definition of the double layer potential operator P require the ex- 
istence of a fundamental solution to div74*V. For divergence form systems, such 
fundamental solutions were constructed by Hofmann and Kim [7] under the hypoth- 
esis that solutions to div AVu = and div A*'Vu = satisfy De Giorgi-Nash local 
Holder estimates. That solutions to divAVu = satisfy such estimates means that 

(1) ess sup,,,,„„„, J"W - < fl-C^)/^ ( / \uA 

\y ~ ^\ \Jb{x;2R) J 

holds whenever m is a weak solution to divAVw = in -B(x; 2i?) C R^'^", for some 
a > 0. It is known that ([T]) is equivalent to the gradient estimate 

(2) \ |Vu|' < (r/i?)"-i+2M /" |y^|2^ < r < i?, 

JB{x,r) Jb{x,R) 

for all weak solutions u to div AVu = in i?(x; R) C R^"*"", for some fi > 0. 

It is known that ([H), or equivalently ([2]), holds for all divergence form systems 
divAVu = where A is real and scalar, and small Loo-perturbations of such {t- 
independence of A is not needed here). Estimates ([1]) and ([2]) also imply the Moser 
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local houndedness estimate 

(3) ess supy,^(,^^) \u{y) \ < (I \u 

\Jb{x;2R) 

whenever divAVu = in B{x.]2R) C R^"*"". We refer to [71 Sec. 2] for further 
explanation of these resuhs. 

At the 8th International Conference on Harmonic Analysis and Partial Differential 
Equations at El Escorial 2008, S. Hofmann formulated as an open problem whether 
(F) as a special case implies the result (S). Our main result in this paper is that this 
is indeed the T) = T) whenever V is defined. More precisely, we prove the 

following. 

Theorem 1.1. Let n,m > 1, and let A = A{x) G Loo(R"; £(C"*(^+"))) be t- 
independent and accretive in the sense that there exists k > such that 

(4) Re/ {A'^{x)r{x)J\x))dx>K I \f{x)\^dx, 

for all f E L2(R"; C'"(i+")) such that curln/n = 0. 

Assume that whenever u is a weak solution to div AWu = in a ball i?(x; 2R), u 
is almost everywhere equal to a continuous function and the Moser local boundedness 

estimate ([3]) holds. Then there exists a fundamental solution T (^t,x) ^ ioc(f^"'^^"! ) 
to divA*V with estimates 

[ \VT^tMs,y)\'dy<iR+\s-t\)-\ 

J\y-x\>R 

for all R > 0, t, s E H and x G R". Moreover 

(5) / {-d,,S'(:^)iO,y),h^iy))dy=ik{BD)h)\{x) 

holds for almost all (t,x) G T{}^'^ and all scalar functions h G L2(R"';C'^). The 
right hand side is defined in Section [3 In particular, we here identify h with a 
normal vector field h G L2(R"; C™'*^^"'"'^^). 

This theorem allows us to transfer known results for the double layer potential 
operator 

Vfh' = Vth' := {ht{BD)h)\, t > 0, 
defined through functional calculus, to the double layer potential operator 

Vfh' = Vth':= [ {-d,^S'(;My),h^iy))dy, t > 0, 

defined classically as an integral operator. The following is a list of such known 
results for Vf, which therefore also hold for = Vf under the hypothesis of 
Theorem ll.il These results for follow by inspection of the proof of [31 Thm. 2.3 
and 2.2], and extends the results for from [H [6]. 
• We have estimates 

snp\\Vth\\l + / \\dtVMltdt+\\N,{Vth)\\l<\\hf, 

t>0 Jo 

for any system with bounded and accretive coefficients A, where the modified 
non-tangential maximal function A'^* is defined in Section [2l In presence of 
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Moser local boundedness estimates of solutions, A^^, can be replaced by the 
usual point wise non-tangential maximal function. 

For any system with bounded and accretive coefficients A, the operators 
converge strongly in L2, that is there exists an L2(R"; C™) bounded operator 
V such that 

lim \\Vth - Vh\\2 = 0, for all h G L2(R"; C"'). 

t—>-o+ 

The map 

{accretive A e Loo(R"; £(0™^'+")))} 3A^V^e /:(L2(R"; C")) 

is a holomorphic map between Banach spaces. In particular, e £(L2(R"; 
depends locally Lipschitz continuously on A G Loo^R"", C{C"^^^~^^^)), and 
therefore invertibility of is stable under small perturbations of A. 
The operator E £(L2(R"; C")) is invertible when A is hermitean, (A^^)* = 
A^^, when A is constant, A{x) = A, and when A is of block form, A^^ = = 

It is also known that is not invertible for many A, even for real and scalar 
(but non-symmetric) coefficients A in the plane, n = 1. A counter example 
was found in [HI Thm 3.2.1] among the coefficients 

1 A;sgn(x) 
-A;sgn(x) 1 



A{x) 



Note that V"^ = V"^ for all these coefficients by Theorem 11.11 It was shown 
in [S] that is not invertible for these coefficients when k = 1. Moreover, 
from [S] and [HI Rem. 5.4] it follows that is invertible for these coefficients 
when k ^ 1, but that the coefficients with k > 1 are disconnected from the 
identity A = I hj the set of coefficients for which is not invertible. 



In the process of proving Theorem 11.11 we also give a new construction of funda- 
mental solutions to divergence form systems. As compared to [7] , this works also in 
dimension 2 and it does not involve any limiting argument but constructs the fun- 
damental solution directly using functional calculus. Extending this construction to 
t-dependent coefficients, we prove the following result. Note that we formulate this 
result in dimension n, not 1 + n. 

Theorem 1.2. Let n > 2, m > 1. Assume that Aq E Loo(R"; >C(R"')) are real 
and scalar coefficients, identified with a matrix acting component- wise on f E C™", 
which are accretive in the sense that there exists k, > such that 

Re(Ao(x)/,/) > K\f\^ for all f E C™",x E R". 

Then there exists e > such that whenever A E Loo(R"; /^(C"*")) is such that 
ess sup^^-^„\A{x)—Aq{x)\ < e, then there exists a fundamental solution to divAV , 

( X — u 

I.e. a function V.^ E Wl^^^{YC; ) such that <1iyA'^W{^ = <^ ^' ' ?n dis- 



tributional sense, with estimates 

(6) f \VTM?dy<R'-, 

J R<\v-x\<2R 



i 7^ k 
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for all R> and x G R". 

The main difference between our construction here and the construction in [7j, is 
that we construct the gradient VF^. instead of the potential r^.. An advantage of 
this approach is that it works in dimension 2 as well. From the gradient estimate 
(E]), we deduce pointwise estimates of Tx in Section [6l This section also contains 
the proof of Theorem 11.21 which builds on the proof of Theorem 11.11 which is in 
Section [5l Sections O [3] and H] contains the details of the construction of the funda- 
mental solution for t-independent coefficients, which uses the Green's formula from 
Definition 13. 11 Half of this identity yields the representation formula for the dou- 
ble layer potential operator. By a duality argument we also derive corresponding 
results for the gradient of the single layer potential operator in Section [7] . 



2. Functional calculus for divergence form equations 



In this section, we explain the method of functional calculus (F) for the L2 Dirich- 
let problem for the equation divAVw = in R^fj"". We assume in this section that 
the coefficients A E L^{lV-~^"'; C{C'^^^~^^'^)) are t-independent and accretive in the 
sense of (jlj). 

Recall from complex analysis the following two relations between harmonic func- 



tions and analytic functions in C = R 



u is harmonic if and only if / = Vu is 



anti- analytic, that is divergence- and curl-free, and (b) u is harmonic if and only if 
there exists an analytic function v with Kev = u. In this section, we generalize this 
result to solutions to div AVu = in R^"*"", following [3l Ej. Following the notation 
from these papers, we shall suppress indices j = 1, . . . , m in this section. 

0, write / = [/x,/ii]* := Vyw]*, where [a, 6]* := ^ 



(a) If divAVw 
Decomposing the matrix A as 



A{x) 



A„ , (x) 



A, 

An 



[X) 

ix) 



we have the conormal derivative di,^u 



A-Aih 



Aj^j^dtu + y4x||V||M, or inversely dtu 
Ax||/||). In terms of /, the equation for u becomes 



dth + diV|, (A||xA:1(/x - Ax„/,|) + Aiiii/ii 



0. 



The condition that / is the conormal gradient of a function u, determined up to 
constants, can be expressed as the curl-free condition 

9,/„ = v„(a:1(/,-a,„/ii 

curl||/|| = 0. 
In vector notation, we equivalently have 



7x" 


+ 





diV|| 






-V|| 








A-l^ 
Au I A , , 



together with the constraint curly /y 



^iiil ■ 
0. Define 



■ All I A , , A I 



0, 



D :-- 





-v„ 



diV|| 




and 



B :-- 



^lix^xx 



A„ 



~^XX^-L|1 

~ ^||X^XX^i|| 
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SO that the equation becomes 
(7) 



dj + DBf = 



together with the constraint R(-D) = {/ G L2 ; curly /y = 0} for each fixed 

t > 0. (Here and below, R(-) and N(-) denote range and null space of an operator.) 
This equation for /, which is an L2(R"; C™(i+")) vector- valued ODE in t, can be 
viewed as a generalized Cauchy-Riemann system. 



Definition 2.1. The conormal gradient of u is the vector field 



ViiM 



where d^^u 



(y4Vu)j^ is the (inward relative Ri',.'^'^) conormal derivative, 
(b) Another Cauchy-Riemann type system related to divAVu = is 

dtV + BDv = 0, 

where D and B have been swapped. Applying D to this equation yields {dt + 
DB){Dv) = 0, so 

/ := Dv = [diV||?;||, -V\\V_l] 
is the conormal gradient of a solution u to div AVu = 0. Looking at /y, we see that 
we should set 

u := —Vj_. 

Then VyU = /y. Moreover 

dtu = -dtv^ = (BDv)^ = (Bf)^ = A-^lih - ^xi|V||u), 
so that di,^u = /x. Thus the equation 
(8) dtV + BDv = 



for V 



implies that u solves divAVu = 0. The vector-valued function v« 



can be viewed set of generalized conjugate functions to u. 

Definition 2.2. A conjugate system for m is a vector field v solving dtV + BDv = 
such that 

u = — fx- 

We now consider the closed and unbounded operators DB and BD in the Hilbert 
space L2 = L2(R'^; C"**-^"*""^). Here D is a non-injective (if n > 2) self-adjoint 
operator with 

R{D) = n and N{D)=n^, 

whereas is a bounded and accretive multiplication operator just like A. Indeed, 
in [3] it was noted that the transform 



A 



A, 
A, 



A, 
A, 



A- 



-A 1 1 ^4 



^l|i^±± ^llll ~ ^l|i^±±^-L|l 



has the following properties. 

(i) If A is accretive, then so is A. 

(ii) If A = B, then B = A. 

(iii) li A = B, then A* = NB*N, where N :-- 
for vectors across R". 



-I 
/ 



is the refiection operator 
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As B is bounded and accretive, we have 

u := sup I arg(5/, /)] < tt/2. 

f&H\{0} 

The operators DB and BD both have spectrum contained in the double sector 

5^_U{0}U5^+, 

where 5"^+ = {A G C \ {0} ; | argA| < u} and S^j- := —5'^+. There are decom- 
positions of L2 into closed complementary (but in general non-orthogonal) spectral 
subspaces associated with these three parts of the spectrum. For DB we have 

L2 = E^L2 © B^'n-^ © E+L2 

and for BD we have ^ ^ 

L2 = £212^71^ (BE+L2. 

Note that for DB we have R(D5) =71 = E^L2 © ^^La and H{DB) = B^^U^, 
whereas for BD we have R{BD) = BH = E2L2 © ^jLg and N(5D) = n^. 
The proof of the fact that the the projections E^ and E^ associated with these 
splittings are bounded uses harmonic analysis from the solution of the Kato square 
root problem. 

Important in this paper are the following intertwining and duality relations. 
Proposition 2.3. We have well-defined isomorphisms 

B : E^L2 E^L2, 
and closed and injective maps with dense domain and range 

D : E^L2 ^ E^L2. 

We also have a duality 

{El,NEi), 

that is the map E^tL2 — ?■ (E'^La)*, mapping g E E^,L2 to the functional E^L2 3 
f I— i- {g, Nf) G C, is an isomorphism. 

Proof. The intertwining by 5 is a consequence of associativity B{DB) = {BD)B, 
the intertwining by D is a consequence of associativity D{BD) = {DB)D, and the 
duality is a consequence of the duality 

{g,N{BD)f) = {g,-{NBN)DNf) = {{-DA*)g, N f). 

□ 

To solve Equation ([7]) for / G "H, we note that DB restricts to an operator in 
E^L2 with spectrum 

a{DB\j,±J C 

Thus e-^^^f is well defined for / G E^La if t > and for / G E^Ls if t < 0. 

The following result was proved in [2] . Here the modified non-tangential maximal 
function of a function / in R^"^" is the function N^:f on R" defined by 

iVJ(x) := sup \t\-^'-''''^/'\\f\\L,(w(t,.)), X G R^ 

±t>0 

where the Whitney regions are W{t,x) := {{s,y) ; Cq^ < s/t < cq, \y — x\ < Ci|i:|}, 
for some fixed constants cq > l,ci > 0. Also, here and below, we write ft{x) := 
f{t,x). 
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Proposition 2.4. Let /q G o^nd define 

f{t,x) := (e-*^^/o)(a:), ±t > 0,x G R^ 

T/ien 

(i) / = [d^j^u, VxuY for a weak solution u to divAVn = in unique up 
to constants, 

(ii) ±(0, oo) 3 t ft E L2 is continuous, with limf_j.o± ft = fo O'^'^d limi_^-|-oo ft = 
in L2 sense, and 

(ill) we have estimates 

II/0II2 ^ sup ll/ill^ ^ // \dtf{t,x)\hdtdx^ [ \N4f)\'dx. 

Conversely, if u is any weak solution to divAVu = m Ri"", with estimate 

||-^*(/)||2 < 00, or sup_|_j>o ll/tib < 00, of the conormal gradient f = [di,^u,'Vxu], 
then there exists fo G E^L2 such that f{t,x) = fo){x) almost everywhere in 

Similar results apply to Equation ([8]). The following result was proved in [2]. 

Proposition 2.5. Let vq G E^L2 and define 

v{t,x) := {e-^^"vo){x), ±t>0,xGR". 

Then 

(i) u := —v^ is a weak solution u to divAVu = in 'R^"', 

(ii) ±(0, 00) 9 t H-> t>( G L2 is continuous, with lim(_^o± t^t = and limt_!.-|-oo Vt = 
in L2 sense, and 

(iii) we have estimates 

\\vo\\l ~ sup \\vt\\l ~ / / \Vu{t, x)\'^tdtdx ^ / \N^,{v)\'^dx. 

±t>0 J JR^+" JR" 

Conversely, if u is any weak solution to divAVu = m R^"^", with estimate 
JJ^i+n \\/u{t, x)\'^tdtdx < 00, then there exists Vq G E^L2 and a constant c G 
such that u{t,x) = —{e~^^^vo)±{x) + c almost everywhere in R^"^"'. 

3. Green's formula on the half space 

Recall that for the Laplace operator, that is the special case A = I and m = 1, 
we have the fundamental solution 

(„_l)/2 



<^{t,x) 



(n-1) 



L-ie + l,fp-'«^, n>2, 



2^ In v^t^ + x^, n = 1, 

with pole (0, 0), where cr„ denotes the area of the unit sphere in R^"*"*^. We note that 

, , 1 

for n > 1. 

In this section, we construct a fundamental solution to more general divergence 
form operators divAV using functional calculus. We assume in Sections [3], H] and [5] 
that the coefficients A G L^(R}'^'^; C{C^^^^'^^)) are t-independent, accretive in the 
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sense of (jlj) and that solutions to divAVn = satisfy the Moser local boundedness 
estimate 1^. 

To explain the definition, we start with the following formal calculation. Assume 
that r = (rj^^^^^^)(t,x))™. is a fundamental solution to divA*V, that is 

Assume that > and that u solves div AVu = in Rj',.^". With appropriate 
estimates of u and F*, Green's formula shows that 



Jr." ^ ' 
where the conormal derivative is d^^^u^ = {A^^Vu'')^. Now let w be a conjugate 
system for u so that = — and dp^u^ = div\\vj. Then by integration by parts, 
we obtain 

u'{to,xo) = - [ (VA.rj;^„)(0,x), W(0,x))rfx, 

where the conormal gradient is V^'T-'* = [diy^^T^^, VnF-'*]*. More generally, it follows 
in this way that if Vq G E^L2, then 



and if Vq G EZL2, then 



{'^A*rl„.o)(0,x),Nvi{x))dx 



(e-*°^^i;o)i(a:o), > 0, 
0, to < 0, 



We now reverse this argument, taking these four formulae as definition. From the 
Moser local boundedness estimate ([3]), it follows that 

\uito,Xo)\ < N^u{x) for \x - Xo\ < Ci\to\/2. 

Thus 

|(e-*«^%)x(xo)|<|tor"/1^o||2, 

uniformly for vq G E^L2 and ±to > 0. Proposition 12.41 and the duality from Propo- 
sition [213] enable us to make the following construction. 

Definition 3.1. For (to,xo) G R^+" and 2 = 1, . . . , m, let Fj.^^^^) = (Fj^^^^))^- be 
the, unique up to constants, weak solution to div/l*VF*^^^ = in Ri"^*^ such that 



for all Vq G E\L2 



For (to, Xq) G RL"*"" and i = 1, . . . , m, let F*^^^ ^.^^ = (F|*^ be the, unique up to 
ution to divy4*VFj, ^ . = in Ri+" such that 



'R" 

for all Vq G E^L2. 
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Some straightforward observations are the following. 

2 

Lemma 3.2. FortQ,a > 0, there is a constant c G such that 

for almost all {t,x) E RL"^". Simlarly, for fixed to < 0, a > 0, there is a constant 
c G C™^ such that 

^ito,xo){t + a,x) = r(t(,_„,^o)(t, x) + c 
for almost all (t, x) G R+^'". 

Furthermore, there are estimates || Vr(tQ^a,Q)(t, ■) ||2 ^ |t — tol""''^- 

Proof. Fix to, a > and consider the functions Va^T*^^^ ^^-^{t—a, ■) = e~'^*~"^'^-^VA*r*^^^ ^ 

and VA*q,„+.,.„)(t, ■) = e-*^^V^*rJ,„+,,,„)(0, ■) in E^^L^, where 5 := NB*N. We 
have 



(VA.rJ^,.„)(t-a,x),iVt;^(x))rfx= / {V A^Tg^^,){0, Ne^'-'^^^''vlix))dx 

R" JR" 

(VA.rj;^^^,^)(0, ■), Ne'^%iix))dx = I (V^.rj;^,^,^)(t, x), Nvi{x))dx 

'R" JR" 

for all Vq G ^j;L2, and therefore ^ A*^\t^^^^){t -a,x) = V A''^ltg+a,xo)(^^ 

The proof for to < is similar. The estimate of || Vr(tQ a;(,)(t, ■)||2 follows from 
Proposition 12.31 and the bound |to|~"^^ of the functional Vq i— )■ (e~*°^^t'o)±(^o)- D 

Note that the translation invariance from Lemma 1X2) enables us to define, for any 
(to,xo) G R^+", a weak solution T(^to,xo)(tyX) to divA*'VT(^to,xo) = in {(t,x) ; t ^ 
to}, so that 

r(to,a:o)(^, 2^) = r(tQ+a,a:o)(^ + a, x). 

We shall prove in the following sections that for appropriate choices of constants, 
this defines a fundamental solution Y [t^^^^^-^it^x^ to divA*V on R^+", that is that 
the traces at t = to coincide except for a Dirac delta distribution at (t,x) = (to,Xo). 
Note that in this paper, except in Section ini we only define the fundamental solution 
on R^'^" modulo constants. 

4. Green's formula on Lipschitz graph domains 

In this section, we improve the estimate || VA*r(j(,,2;Q)(t, ■) ||2 ^ |t — to|~"^^ away 
from Xo, and prove the following. 

Proposition 4.1. We have for R> and t ^ to the estimate 

[ \VA'T^t,,^,){t,x)\^dx <iR+\t- to|)-^ 

J \x—xo\>R 

To prove this, we consider the graph 

S = {(7(x),x) ; X G R"} 
of a Lipschitz function 7 : R" — R. We assume 7(xo) = and write 

Vi|7Glvoo(R";R"). 
Recall the following consequence of the chain rule. 
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Proposition 4.2. Let Vt C R^"^" he an open set. Then u is a weak solution to 
div AVu = in {{t + 'j{x),x) ; (t, x) G Q} if and only if 

Ua{t, x) := u{t + 7(x), x) 

IS a weak solution to divA^Vu^ = inVL. Here 



1 -a' 
I 



1 
-a I 



has estimates \\A„\\oo ^ (1 + IkllDll^lU and ka, > ^^/(l + ||cr||^), where a* 
denotes the transpose of the column vector a. 

Proposition 4.3. Fix (to^xo) E R^"^" and consider a Lipschitz graph S as above 
such that 7(x) > and 7(xo) = 0. Define T = T(^tg^xQ) for coefficients A, and define 
To- = ^(to,xo) for coefficients A„, as in Definition Then there is a constant 

c G C™^ such that 

for all t <to, X e R"-. 



r{t, x) = r„{t - 7(x), x) + c 



Proof, (i) The function r(t,x) is uniquely, up to constants, determined by tlie prop- 
erty tliat 



(9) 



R" 



{VA'T^\0,x),Nvl{x))dx = (e^*»^^t;o)i(xo). 



for all vq G E\L2. By tlie intertwining B : E\L2 — )■ E\L2 from Proposition 12.3^ we 
can write 

for a weak solution u to div AVu = in RV^". Then dHD reads 



(10) 



R" 



VF'(0,x) 



-Ai\ 
Alt 



VM^'(0,a;) ) dx = dtu\to,xo). 



vri^(o,x) 







V<(0,x) dx = dtul{to,xo). 



(ii) Now replace A, T and u hj Acr, To- and 

Ua{t, x) := u{t + 7(x), x). 

Then Mo- is a weak solution to divA^Vw^ = for t > —7(0;), and in particular 
Va,m^(0, ■) G ^X^2. As in (Uni), we have 



R" \ L " ^"'-'J 

Here V„ is a weak solution to div(A^*)*Vr^^' = for t < to- Since {A„)* = (A*)^, 
the function 

T(t,x) := T„(t - j{x),x) 

is a weak solution to div(A''*)*Vr''^ = for t < + 7(3;)- Changing variables in 
Til) , we get 

VF"(7(x), x), A^'*^(-l, a(x))VM*^(7(x), x))rfx = ^^^'(to, xq). 



(12) 
where 
(13) 



R" 



Ai'=(zyo,zy) := A^'* 



'^^0 0" 




z/* " 


z/ 


+ 


-z/o_ 
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This uses the chain rule 
and the calculation 



'l 


o" 


t - 


a 


I 









1 0' 

a{x) I 







VM'=(t + 7(x),x), 





'l 


o" 




a 





A^-^(-l,a(x)). 



(iii) We now apply Stokes' theorem to the 1-form n-form 

on {{t,x) ; < t < 7(x)}. Using ( fT3l) and the product rule shows that its exterior 
derivative is 

(yv^\t,x),K^\duW^;)Wu\t,x)) = {dw{A^''ywv^\dtu^) + {{A^^ywv\dtWu^) 

+ {W^^dtV^' - dtW{r^\ (A^''=Vu'=)||) + (VF\ [diV||(A^''=Vu^)||, -9t(A^''=Vn'^)||]*) 
= + {{A^^)*Vr\dtVu^) + + (VF\ -^^(A^'^Vm^)) = 0. 
Thus, applying Stokes' theorem to ffT^ gives 



dtu\to,xo) = I {VP\0,x),A^\-l,0)Vu\0,x)]dx 

vr%o,x 

Comparing with (ITOll proves the proposition, by uniqueness of T 



R" 



-A^A ■ 

Ai' 



Vu''{0,x) dx 



□ 



Proof of Proposition \4.1\ Consider T(^tg^xo)- By Lemma |3.2[ we may assume that 
t = 0. Fix R > to > and apply Proposition 14.31 with 



7(x) 



Then || Ao-||oo ~ ll^lloo and ^ ka and r(0, x) = Tfj{—R, x) + c for > R. Thus 
the estimate from Lemma [3.21 gives 



1 kl> 


\x\ 


[R, 


\x\ 



/ \VA*T{0,x)\'^dx = / \VA*Ta{-R,x)\'^dx 

J\x—xo\>R J\x—X(,\>R 



<||v^;r.(-i?,-)f <|i? + ^o|-^ 

This proves the estimate for > t- The proof for to < ^ is similar, using the analogue 
of Proposition 14.31 for to < t. □ 

5. Fundamental solution for t-iNDEPENDENT coefficients 

In this section, we complete the proof of Theorem 11.11 Fix xq G R" and i = 
1, . . . , m, and define the vector field 

f (t,a;):=V^.rg,^)(t,a:), t^O. 

As in Section [21 we suppress the index j. 
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Proposition 5.1. For R > 0, we have the estimate 



\f{t,x)\'dtdx < R' 

R<\{t,x)-{0,xa)\<2R 

In particular, for 1 < p < {n + ^)/n, we have f G L^°^ (R^'^'^) and 

ff \f{t,x)\^dtdx<R'-<^-^\ 

J J \{t,x)~{0,xo)\<R 

Proof. From Proposition 14.11 we obtain tlie estimate 

r^^ dt 

\f(t,x)\^dtdx< / — ^ ^<^'~"- 

' B<\it,x)-io,xo)\<2R Jo (max(i? - t, 0) + t)" 

Holder's inequality then gives the Lp-estimate after summing a geometric series. □ 
Proposition 5.2. We have that f G L^°'^ (R}'^"') and, in R^"*"" distributional sense, 

('5(o,xo),0), j = i, 
,0, j^t, 

and curl||/|| = 0. 



m + DB)fy 



Proof. That curl||/jj' = is clear from the construction of /. To compute {dt + DB)f, 
we fix a test function G Co^(Ri+'^; C'"(i+")) and define 



I{xo) := I / (/(t, x), -N{dt + B{x)D)(j){t, x))dtdx. 

/Rl + n 



For e > 0, let 



Since 



J,(xo) := / / (/(t, x), -N{dt + B{x)D)<p{t, x))dtdx. 
'|i|>. 



|J(xo)-/.(xo)| < // \f{t,x)\dtdx 

J J \t\<€,\x\<C 

for some C < oo depending on </>, we have that le ^ I uniformly in xq as e — )■ 0. By 
Definition 13. 11 we have 

roc ^ p—e 

/,(xo) = - / (e*^^E^(5D + dt)<i)trAx,)dt + / {e'''''E\{BD + 9i)0i)l(xo)rft 

J e J — oo 

/•oo /* — e 

>/ — oo 

Therefore G L2 and 

in L2 as e — )■ 0. We have here used that = for v G N{BD). 

We note that /e are continuous functions, since A satisfies property (M), and 
converge uniformly to /. Thus / is continuous, and it suffices to prove 



|/(xo) -0(O,xo)i|^c/x = 

K 



LAYER POTENTIALS BEYOND SIOS 



15 



for an arbitrary compact set K. But this is clear since — (0o)x ^"^^ — )■ / in 
L2{K). This proves the proposition. □ 

Proof of Theorem \l.l[ Fix (tg, Xq) G R^^*^, define as in SectionOthe function T(^tQ,xQ){ii x) 
for t 7^ to- Without foss of generahty, we can assume that Iq = 0. It foUows from 
Propositions 14. 1^ 15.11 and 15.21 that r((g x) is a fundamental solution with the 

stated bounds, using the correspondence between divA*V and dt + DB from Sec- 
tion H 

By Definition 13. we have for all Vq G E\L2 and almost all (t, x) G 'R}^'^ the 
identity 

/ {VA'T\,^4Q,y),Nvo{y))dy = (e-*^%)l(x), 

where r = Now let vq = E\h = E^[h,OY for some scalar h G L2(R'^; C™), 

or equivalently normal vector field [h, 0]* G L2(R"'; C™'^^"'""-'). We then obtain 

/ {VA*ri,^,^{0,y),NEXh{y))dy = {e-'^^E\hUx). 

JR" 

Using the duality from Proposition 12.31 the left hand side is 

/ {E-^,VA*T\,^.){^,y),Nh{y))dy= [ {V A'T\,^40,y), Nh{y))dy 

= - [ {d.,S\t,40,y),h{y))dy, 
whereas the right hand side is {ht{BD)h)\{x). This proves the theorem. □ 

6. Fundamental solution for t-DEPENDENT coefficients 

In this section, we prove Theorem 11.21 and show some further estimates of the 
constructed fundamental solutions. We assume throughout this section that n>2 
and m > 1, and that Aq and A are as in the hypothesis of Theorem 11.21 where we 
choose e > small enough so that the De Giorgi-Nash local Holder estimates ([T]), 
or equivalently ([2]), hold for A- and for y4*-solutions, and that A is accretive. Note 
that in this section we allow Aq and A to depend on all n variables. As in the proof 
of Theorem II. H we write F* = (F-^*)j and suppress the index j, and sometimes also 
i. 

Proof of Theorem \1.2 . (i) Define, in R^"*""-, t-independent coefficients 

i(t,x)[/„/„r := [/x,A(x)/„], 

so that A{t, x) = A{x). Our aim is to construct a fundamental solution for A on R"' 
from the already constructed fundamental solution for A on R^"*"*^, by integrating 
away the auxiliary variable t. We assume that e > is small enough so that A is 
accretive and that De Giorgi-Nash local Holder estimates ([1]), or equivalently ([2]), 
hold for A- and for yl*-solutions. 

In particular this means that the hypothesis of Theorem 11.11 is satisfied for A*, 
giving a fundamental solution F(o,a;o) ^ ioc(-f^" i C™) with pole at (0, xq) to divAV 
in R^"'""' with estimates 

[ \Vf^o,xo){t,x)\'dx<{R+\t\)-\ 

J \x—xq\>R 
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for all /2 > 0, t e R and Xq G R''. 
(ii) Assume first that n > 3. Define 

gl{x) := Vi,J;o,.o)(^'^)' 
so that \\gt\\L2{R<\x\<2R)dt < (^4)"/^ ~ Thus 

/■oo 

:= / i9lix)),dt 



converges in L2{R < \x\ < 2R), and we have ||fi'||L2(R<|a;|<2/?) ^ "''^ so that 

g e L'°=(R"; C™"). It suffices to show that (divA^*)-'' = \ " and curl^^ = 

[0, 3 ^ I, 

in R"-distributional sense. The latter is clear from the definition of g. To prove the 
former, let G C^ilV; C™). Let r] e C^(R) be such that = 1 for |t| < 1 and 
?7 = for > 2, and let ?7t(^) := vi't/T). Consider the integral 



It:= / / {gl{x),A*{x)VtA<l>i^)VTmdtdx = -<j)\xo). 

Then 

h = [f {{gl{x))^,(P{x))dtVT{t)dtdx- [ [ {{gl{x)\,A*{x)V(P{x)){l - r]T{t))dtdx 



{g\ A*V(P)dx =: IIt - Uh + j {g\ A*V4>)dx. 

The estimates \\gt\\2 ^ t~"'^'^ proves that IIt and 1 1 It — j- as T — )■ oo, so that 
J{g\A*Vc(>)dx = -^'{xo). 

(in) Now let n = 2. We claim that in this case 



/ \\9t\\L2{\x\<R)dt < oo. 

J\t\>R 



sup 

R>0 J\t\>R 

From this claim, it will follow that 

,, ,, ^ /-^^ dt P I, „ , ^ , 

\\g\\L2{R<\x\<2R) ^ / D , + + / \\9t\\L2{\x\<2R)dt <1 

Jo -n^ + 1 J2ij 

and IIt and 7//^ — !■ as T — oo as in (ii). To prove the claim, we apply 
the estimate ([2]) to the solution divAVV^^ = in {(t,x) ; max(|a;|, \t — T\) < 
\T\/2} D {{t,x) ; max(|x|, \t-T\) < R} for |r| > 2R. We obtain 

fT+R 

II _ II Jj. ^ . /Till _ll ^ . T-> /rn\l/2+u\\ ,J\ 

2{\x\,\t-T\<\T\/2) 



/ \\9t\\L2{\x\<R)dt < y/R\\g\\L2(\x\,\t-T\<R) < V^(-R/^)^^^^^II5'I|l2( 
Jt-r 

< ^/R{R/\T\Y"+^ ( [ t~'dt) ' < (i?/|T|)i+^ 

V^|t-T|<|T|/2 / 



From this it follows that 

oo 



/ \\9t\\L2i\x\<R)dt <y2{R/kRy+^ <i. 

J\t\>R i^^^ 

This completes the proof of the theorem. □ 
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Proposition 6.1. Under the hypothesis of Theorem and suitable choices of 
integration constants, the following holds. 

(i) The gradient of the fundamental solution to divAV has estimates 



[ \VTM\'dy<r^-'+''^\z- 

J B{z,r) 



X 



for 2r < I2; — x| and some /i > 0. 
(ii) If n > 3, then the fundamental solution to divAV has pointwise estimates 

|r.(y)|<|y-x|2-", y^x, 

and Holder estimates 

\T^iy') - TM\ < \y - ^r-", W -y\<\y- x\/2. 

(ill) If n = 2, then the fundamental solution to divAV has pointwise estimates 
|rx(z/)| < 1 + |ln|i/-x| |, y^x, 
and Holder estimates 

\V,{y') - VM\ < (^)" (1 + I In |y - x| I), \y' - y\ < \y - x\/2. 

Proof, (i) This follows by from and 
(ii) For R> 0, consider the mean values 

{A — i)an-in jR<\y-x\<2R 

We obtain from Poincare's inequality, with means over the inner/outer halves of the 
annuli, and (EI) that 



{VM-AR)dy 

R<\y~x\<iR 



1 



{VM-MR)dy 

R<\y~x\<AR 



< 



< ^ / - Ar\ + \TM - A2R\)dy 

^ jR<\y~x\<4:R 

1 /2 

f ^ / ilUy) -Ar\' + \Tx{y) - A,R\')dy) 

\^ J R<\y~x\<AR J 

< fi?'^" / \VTxiy)\'dyY' <R'-\ 

\ JR<\y~x\<AR J 



If n > 3, we obtain the estimate 

(14) <(2'^)2-", 

for all J, G Z with j > k. In particular limj^oo exists. Choosing the constant 
of integration, we assume that this limit is zero. This gives 



Ar\ < i?^-", for all > 
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and again by Poincare's inequality and that 

[ l^Ml'dyV' <\Ar\+(^ [ \TM-AR\'dyY' 

\^ jR<\y-x\<2R / jR<\y-x\<2R / 

< \An\ + (r'-- [ \VTM\'dy) ' < R'^''. 

\ Jr<\v-x\<2R J 



' R<\y-x\<2R 

Using the Moser local boundedness estimate ([3]) and the De Giorgi-Nash local Holder 
estimate ([T]), this proves the estimates (ii). 
(iii) If 77, = 1, the equation ffT^ becomes 

\A2j - A2k\<i - k. 

Choosing the constant of integration so that \Ai\ < 1, this gives 

\An\<l + \\nR\. 

The point wise estimates (iii) then follows as in (ii). □ 

7. The gradient of the single layer potential operator 
We end this paper by deriving results for the single layer potential operator 

Sth\x)= [ Tl {t,x)h^{y)dy, 

where F here denotes the fundamental solution for divAV. Recall that the Neumann 
problem, with boundary datum ip, is solved through the ansatz u(t,x) := Sth{x), 
where the auxiliary boundary function h solves the integral equation 

lim d^^Sth = ip. 

We prove the following result for the single layer potential operator, analogous to 
Theorem 11.11 for the double layer potential operator. 

Theorem 7.1. Assume the hypothesis of Theorem with A replaced by A*, so 
that r now denotes the fundamental solution for divAV. Then 

(15) Va [ Tl^^{t,x)y{y)dy = {e-'^''E\h)\x) 

holds for almost all (t, x) E iV^"' and all scalar functions h E L2(R"'; C™). We here 
identify h with a normal vector field h E L2(R"; C"^'-"'^"'""^) on the right hand side. 

This theorem allows us to transfer known results for the conormal gradient of the 
single layer potential operator 

VA§fh' = VASth' := {e-'^^E\h)\ t > 0, 

defined through functional calculus, to the conormal gradient of the single layer 
potential operator 

VASfh' = VASth' := / Tl Jt, ■W{y)dy, t > 0, 

defined classically as an integral operator. The following is a list of such known 
results for VA<Sf which extends the results for VyiiS/^ from [H [6] . 
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• We have estimates 

sup||VA5i/i||2+ / \\dt^ j,sMlidt^\\N,{y Asml < \\hf, 

i>0 Jo 

for any system with bounded and accretive coefficients A. 

• For any system with bounded and accretive coefficients A, the operators 
VyiiSj converge strongly in L2, that is there exists an L2(R"'; C™) bounded 
operator Va^ such that 

lim WVASth - VAShh = 0, for all h e L^IW; C"). 

f-i>0+ 

• The map 

{accretive A G L^iK"; £(C™(i+")))} 3 A^ V aS^ G C^L^) 

is a holomorphic map between Banach spaces. In particular, Va^/^ G C{L2) 
depends locally Lipschitz continuously on A G Loo(R"; >^^(C'"^^"^"'^)), and 
therefore invertibility of d^^Sf is stable under small L^o perturbations of A. 

• The operator d^^Sf^ G £(L2(R"; C™")) is invertible when A is hermitean, 
{A^^y = A^^, when A is constant, a4(x) = A, and when A is of block form, 
y4Y|| = = A^j^. The counter example to invertibility mentioned in the 
introduction applies also to du^St. 

Proof. In the classical case of integral operators, the conormal derivative of the single 
layer potential is dual to the double layer potential operator. Similarly, the proof of 
Theorem 17.11 is by duality. We note from Definition 13.11 that 

for all t > and Vq G L2. Integrate this equation against a scalar/normal vector 
field h G L2(R"; C") to obtain 

Since A*D = NB*ND = N{-B*D)N, we obtain 

Since vq G L2 is arbitrary, this proves ( |T5|) . □ 
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